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ABSTRACT 

With their long mean free paths and efficient heating of the intergalactic medium (IGM), 
X-rays could have a dramatic impact on the thermal and ionization history of the Universe. 
Here we run several semi-numeric simulations of the Dark Ages and the Epoch of Reioniza- 
tion (EoR), including both X-rays and UV radiation fields, attempting to provide an intuitive 
framework for interpreting upcoming observations. We explore the impact of X-rays on var- 
ious signals: (i) Reionization history: including X-rays results in an earlier, slightly more 
extended EoR. However, efficient thermal feedback from X-ray heating could yield an ex- 
tended epoch in which the Universe was « 10% ionized, (ii) Reionization morphology: a 
sizable (~10%) contribution of X-rays to reionization results in a more uniform morphol- 
ogy, though the impact is modest when compared at the same global neutral fraction, xhi- 
Specifically, X-rays produce a dearth of fully neutral regions and a suppression of small- 
scale (k > O.lMpc -1 ) ionization power by a factor of < 2. However, these changes in mor- 
phology cannot be countered by increasing the bias of the ionizing sources, making them 
a robust indicator of an X-ray contribution to the EoR. (iii) The kinetic Sunyaev-ZeVdovich 
(kSZ) effect: at a fixed reionization history, X-rays decrease the kSZ power at / = 3000 by 
~ 0.5/iK 2 . Our extreme model in which X-rays entirely drive reionization is the only one 
which is marginally consistent with the recent upper limits on this signal from the South Pole 
Telescope, assuming no thermal Sunyaev-Zel'dovich (tSZ) - dusty galaxy cross-correlation. 
Since this extreme model is unlikely, we conclude that there should be a sizable tSZ-dusty 
galaxy cross-correlation, (iv) The redshifted 21cm signal: the impact of X-rays on the 21cm 
power spectrum during the advanced stages of reionization (xri < 0.8) is modest, except 
in extreme, X-ray dominated models. The largest impact of X-rays is to govern the timing 
and duration of IGM heating. In fact, unless thermal feedback is efficient, the epoch of X-ray 
heating likely overlaps with the beginning of reionization (xri > 0.9). This results in a 21cm 
power spectrum which is ~ 10-100 times higher than obtained from naive estimates ignoring 
this overlap. On the other hand, if thermal feedback is efficient, the resulting extended epoch 
between X-ray heating and reionization could provide a clean probe of the matter power spec- 
trum in emission, at redshifts more accessible than the Dark Ages. 

Key words: cosmology: theory - intergalactic medium - early Universe - reionization - dark 
ages - X-rays: galaxies 



1 INTRODUCTION 

Cosmic reionization is the last global baryonic phase change, and 
acts as powerful probe of early Universe physics and the first gen- 
erations of galaxies. Although the details remain unclear, it is gen- 
erally accepted that reionization is mostly driven by stellar sources 
inside dwarf galaxies that accrete gas via molecular or atomic hy- 
drogen line cooling. The relatively soft ultraviolet (UV) spectra of 
such stellar sources results in a reionization which is "inside-out" 
on large scales: ionized bubbles grow around the first, highly clus- 
tered sources, with their eventual coalescence signaling the com- 
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pletion of reionization (e.g. Furlanetto et al. 2004, McQuinn et al. 
2007 ; Trac & Cen 2007 ; Zahn et al. 201 1). Such UV-driven reion- 
ization scenarios result in a very inhomogeneous ionization field, 
where the large-scale clustering of galaxies determines the ioniza- 
tion morphology. 

On the other hand, reionization would be noticeably different 
if it were driven by higher energy, X-ray photons. Firstly, reioniza- 
tion morphology would be more uniform, due to the relatively large 
mean free path for X-rays (e.g. McQuinn 2012| ): 

Ax " 2 °^K3^v) 26 ( i ^ £ )" 2cMpC ' (1) 
Here xhi is the mean neutral fraction of the intergalactic medium 
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(IGM) and Ex is the photon energy. Secondly, most of the X- 
ray energy gets deposited as heat when xm surpasses a few per- 
cent(e.g.[ Shull & va n Steenberg|1985[|Furlanetto & Stoever|2010[ 
Vald es et al.| |2010), quickly raising the Jeans mass in the IGM. 
Therefore, if X-rays were abundant in the early Universe, thermal 
feedback might delay the end stages of reionization (e.g. |Ricotti &| 
Ostriker 2004, Kuhlen & Madau 2005 ). Such a scenario could ac- 
commodate a relatively large optical depth to Thompson scattering 
( Ko matsu et al.|201l] ) and still allow stellar UV emission to com- 
plete reionization at z ~ 6-7, consistent with recent observations 



(e.g.|Bolton etal.|201ll 


Dijkstra et al.|201 l||Pentericci et al.|201 1| 


ISchroeder et al.|2012j>| 





Several observations indirectly limit the allowed contribution 
of X-rays to reionization. Firstly, the number density of bright ac- 
tive galactic nuclei (AGN) decreases sharply towards high redshifts 
(e.g. |Madau et al.|1999l|Fan et al.|200l) . Although faint AGN (so- 
called miniquasars) could still contribute (e.g. |Madau et al.| 2004 ; 
|Volonteri & Gn edin 2009 ), the parameter space for such models 
is decreasing with more sensitive surveys probing further down 
the AGN luminosity function (e.g. |Willott et al.|[20 10). Attention 
has also focused on other sources of X-rays, such as high-mass X- 
ray binaries (HMXBs), which dominate the X-ray contribution of 
star-forming galaxies (e.g. Mi neo et al.|2 012). The strongest con- 
straints on X-ray reionization comes from the present-day unre- 
solved X-ray background (XRB). The XRB is an indirect probe, 
since the observed photons originating from high-redshifts are en- 
ergetic enough to have mean free paths greater than the Hubble 
length. Therefore the observed photons are not the ones interacting 
with the IGM. Nevertheless, estimates of the unresolved XRB from 
the Chandra deep fields (CDFs; |Hickox & Markevitch|2007] > rule 
out reionization models driven by hard X-ray spectra (with photon 
, Dijkstra et al.|2004[|McQuinn|2012| >, pro 
ral break towards high-energies (Ex > 10 



energy indices a < 1; 
vided there is no spectraT1 
keV). 

Despite the limited parameter space, there has recently been 
a resurgence of interest in X-ray reionization. This has been moti- 
vated by several factors. 

(i) The observed population of high-redshift galaxies (z ~ 6- 
10) seem unable to maintain an ionized Universe, baring an evo- 
lution in galaxy properties (for example, an increase with redshift 
in the fraction of ionizing photons capable of making it our of the 
galaxy; e.g. |Mitra et aLpUH) |Kuhlen & Faucher-Giguere|[20T2l 



IFontanot et al.|2012| ). 

(ii) The observed anti-correlation of HMXBs with metallicity 
(e.g. |Crowther et al.|2010l|Kaaret et al.|201 1| ), as well as theoretical 
studies proposing a high binary fraction among the first stars ([Turk] 
|et al.|2 009 ; Stacy et al. 2 0T0]>, motivate a high er contribution from 
HMXBs in early galaxies ( jMirabel et al.|20 1 lp] 

(iii) The recent upper limit on the kinetic Sunyaev- Zel'dovich 
(kSZ) signal from the South Pole Telescope (SPTf]by [Reichardt 



1 However, similar extended reionization scenarios could result without 
X-rays: for example from strong, self -regulating feedback on small galax 



ies (e.g.|Springel & Hernquist|2003[|Ahn et al.|2012[|Pawlik et al.|2012| 



and/or from reionization "stalling" when the HII regions grow to surpass the 
attenuation lengt h of ionizing UV photons jFurlanetto & Mesi nger 2009 
Crociani et al. 2011 Alvarez & Abel 2012). 



z Note however that there is no evidence of such an evolution of the X- 



ray luminosity to star formation rate (SFR) rat io, up to z ~ 6 (|Cowie et al7 
20121. This z ~ 6 constraint is however strongly dependent on the assumed 

SFR and X-ray spectral index. 

3 http://pole.uchicago.edu/ 



|et al.| ( |2012| > is inconsistent with stellar-driven UV reionization 
models ( Mesinger et al. 2012), unless there is a non-negligible cor- 
relation in the thermal SZ - cosmic infrared background (CIB) sig- 
nal. This measurement is a probe of the ionization structure on 
~ 20 Mpc scales; X-rays with long mean free paths are a nat- 
ural candidate for erasing such structure ( Visbal & Loeb |2011| 
IMesinger et al.|2012| >. 

(iv) |Cappelluti et al.| ( |2012| ) recently detected a correlation be- 
tween the CIB and the soft XRB. The observed cross-power is best 
explained if the sources are at high redshifts (z > 7). 

In this work, we do not attempt to predict whether or not X- 
rays are important to reioni zation (see, e.g. |Dijkstra et al.||2004[ 
Srbin ovsky & Wyithe|2007l |Volonteri & Gnedin|2009| |McQ uinn 



2012) . Instead we explore various observational signatures of X- 
rays in the reionization and pre-ionization epochs. We attempt to 
provide an intuitive framework for interpreting upcoming observa- 
tions. Given the many astrophysical uncertainties, we adopt a para- 
metric approach, looking for general trends and physical insight. 

We model X-rays in the early Universe using a public code, 
21CMFAST. This is a "semi-numerical" simulation, using well- 
tested and efficient algorithms to compute various cosmic fields. 
Due to their long mean free paths, very large boxes (several hun- 
dred comoving Mpc) are required in order to properly model X-rays 
(as well as the photons responsible for pumping the 21cm line; see 
below). Computing these radiation fields for the dominant galaxy 
population in the early Universe is computationally impractical for 
conventional radiative transfer (RT) algorithms. On the other hand, 
analytic approaches do not preserve the spatial structure we are 
investigating. Therefore semi-numerical simulations are uniquely 
suited for this task. 

This paper is organized as follows. In ^2] we present our sim- 
ulations. In ^3]we discuss our results, including the impact of early 
X-rays on the evolution of the global reionization signal (j ]3.1| ), 
reionization morphology ( §3.2| ), the kSZ signal ( §3.3| >, and the red- 
shifted 21cm signal ( §3.4| ). Finally, in Q we summarize our con- 
clusions. 

Except for fluxes, we quote all quantities in comoving units. 
We adopt the background cosmological parameters (Qa, Qm, 
Q b , n, cr 8 , H ) = (0.73, 0.27, 0.046, 0.96, 0.82, 70 km s" 1 
Mpc -1 ), consistent with the seven-year results of the WMAP satel- 
lite (WMAP7; [Komatsu et al.|20TT) . 



2 SIMULATIONS 

We use a parallelized version of the publicly available 21CMFAST 
codeQ This code uses perturbation theory (PT) to generate the den- 
sity and velocity fields, and PT plus the excursion- set formalism of 
Furlanetto et al. (2004) to generate the ionization fields. We sum- 
marize these calculations below. For further details on these algo- 
rithms, see [ Mesinger & Furlanetto| ( [2007| ), [Zahn et al.|p0TT] ) and 
Mesi nger et al.| ( |2011| >! 



2.1 Ionization and heating by X-ray photons 

In the early stages of reionization, there is a small HII region local 
to each source (e.g. Cen 2006), whereas the X-rays are able to pen- 
etrate deep into the mostly-neutral IGM. Outside of the local HII 
regions, we compute the temperature, Tk, and ionized fraction, x e , 



http://homepage.sns.it/mesinger/Sim.html 
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at a given redshift, z, and position in our simulation box, x, accord- 
ing to: 



(ix e (x, z) _ dt 



dz 



c?T K (x, z) 
dz 



dz 



[r, 



a A Cx e n b fu] 



(2) 



dt 



3fcs(l + X e ) C/Z 



+ 



2Tk <iri6 Tk dx e 
3rib dz 1 + x e dz 



(3) 



where rib = ^b,o(l + ^) 3 [1 + £ n i(x, z)] is the total (H + He) 
baryonic number density at (x, z), e p (x, z) is the heating rate 
per baryon for process p in erg s _1 , Tion is the ionization rate 
per baryon, q>a is the case-A recombination coefficient, C = 
(n 2 ) I (n) 2 is the clumping factor on the scale of the simulation 
cell (here taken to be C = 2), ks is the Boltzmann constant, fn is 
the hydrogen number fraction. 

For X-ray heating and ionization, we have: 



ex(x, z) 



and 



-/ 

J \ 



Max[^o,^ T = i] 



hv 



ix, z 



I 



Max[z^o ,v T =l] 



(4) 



dV ~Tw ^ fi X i°iFi ) ( 5 ) 



"( 



hv - EX 



rth\ n,HI 



E th 

^HI 



/"ion, Hel ^ ./ioi^Hell 



77-th 
^Hel 



^Hell 



+ 1 



where the frequency integral includes the contribution from pho- 
tons within a mean free path of (x, z J] and we perform a sum over 
species, i = HI, Hel, or Hell, in which fa is the number fraction, 
Xi is the cell's species ionization fraction [which for HI and Hel is 
(1 — x e ), and for Hell is x e ], o% the ionization cross-section, and 
Ef 1 is the ionization threshold energy of species i. Furthermore, 
/heat [hv - Ef 1 , x e (x, z)] and / ion ,j [hv - Ef 1 , x e (x, z) , j] are the 
fraction of the primary electron's energy going into heat and sec- 
ondary ionizations of species j respectively (taken from Furlanetto 
|& Stoever|20T0] ). 

In the above equations, the angle- averaged specific intensity, 
J(z/, x, z), (in erg s _1 Hz -1 cm -2 sr _1 ) can be computed integrat- 
ing along the light-cone: 

jf ^ (i+^) 3 r a > cdt <f* 

and the comoving specific emissivity is evaluated at v e — v{l + 

z')/{l + z): 



€ h v(v e ,x, z) =ah 



Lim p \v 

Pcrit,0^b/*(1 + ^nl) 



df coil 



dt 



(V) 



5 Note that this spatially-averaged, frequency-dependent opacity is com- 
puted including the total ionized fraction from both X-rays and UV photons 
(seelMesinger et al.|2011l. 



where iVx is the number of X-ray photons per stellar baryon, fim p 
is the mean baryon mass, p C rit,o is the current critical density, 
is fraction of baryons converted into stars (we take /* = 0.1), S n \ 
is the mean non-linear overdensity and / co ii is the fraction of mat- 
ter collapsed into halos with virial temperatures greater than T v i r , 
which is computed according to the hybrid prescription of Barkana 
& Loeb ('2008). The quantity in the brackets is the comoving star 
formation rate density (SFRD). The specific emissivity is chosen to 
have a spectral (energy) index of a = 1.5 both in order to facili- 
tate comparisons with previous works (e.g. Pritchard & Furlanetto 
12007) [Santos et al.|2008||Baek et al.|20101|Mesinger et al.|2011| >, 
and also since a significant contribution from harder spectra is ruled 
out by observations of the XRB (McQuinn 2012) unless there is a 
break at high energies (Ex > 10 keV). 



2.2 Ionization by UV photons 

The contribution of UV photons to reionization is computed as a 
post-processing step, following the X-rays. We use the excursion- 
set prescr iption of jFurlanetto et al.| ( |2004| >, slightly modified ac- 
cording to Zah n et al.|p011| ) (our 'FFRT' scheme). This algorithm 
determines whether a region is ionized by comparing the (inte- 
grated) number of ionizing photons to the number of neutral atoms. 
Specifically, a simulation cell at coordinate x is flagged as ionized 
if 



Cuv/coll(x, Z,R)>1- X e (x, Z, R) 



(8) 



where Cuv is an ionizing efficiency parameter (described below), 
and as above, / co ii is the fraction of mass residing in dark matter 
halos inside a sphere of radius R and mass M 4/3ttR 3 p, where 
p — p[i -f S n \]. In this work we include on the RHS the mean 
fraction of neutral hydrogen remaining within R: 1 — x e , computed 
according to eq. ([2]). 

As for the X-rays, this tabulation includes all halos above a 
set virial temperature threshold, T v i r . We do not explicitly compute 
individual halo locations (as in Mesinger & Furlanetto 2007), but 
instead apply the conditional excursion-set formalism directly on 
the non-linear density field to compute / co ii (Zahn et aL] 201 1). The 
ionization field is also computed following the excursion-set ap- 
proach (e.g., |Bond et al.|1991|[Lacey & Cole|1993l|Furlanetto et al.| 
2004). We decrease the filter scale R, starting from some maximum 
value i?max- A physical choice for i? max corresponds to the mean 
free path of ionizing photons, R m f P (Furlanetto & Oh 2005 ; Choud- 



|hury et al.|2QQ9]| Alvarez & Abel|2012) . We take i? ma x 40 Mpc, 
consistent with recent observations (e.g. Songaila & Cowie 2010) 
and theoretical estimates (e.g. McQuinn et al.|201 1\ . If at any R the 
criterion in equation (8} is met, this cell is flagged as ionized. 

Note that this formalism assumes that UV photons, with their 
short mean free paths, carve out fully-ionized HII regions in the 
IGM (e.g. Zahn et al. 20lTJ. If we denote the covering factor of 
these HII regions as Qhii, then the total ionized fraction (includ- 
ing the partial ionizations from X-rays in eq.[2]), can be written as 
Xi ~ Qhii + (1 — Qmi)x e . Similarly, we denote the mean neutral 
fraction as xhi = 1 — x%. 



2.3 The runs 

Our simulation boxes are 750 Mpc on a side, with a resolution of 
500 3 . We integrate the evolution from z — 40 down to z — 5.6, in- 
cluding the end of the Dark Ages and reionization. Here we discuss 
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the astrophysical parameters we vary, and motivate some fiducial 
choices. 

We parameterize the ionizing efficiency of UV photons as fol- 
lows: 



Cuv 



iV 7 
4400 



fesc 

oT 



o.i 



1.5 



1 + n T , 



(9) 



where / esc is the fraction of UV ionizing photons that escape into 
the IGM, iV 7 is the number of ionizing photons per stellar baryon, 
and n rec is the mean number of recombinations per baryon (a fac- 
tor which appears since Cuv is defined after integrating the ion- 
ization rate over cosmic time). A value of /uv ~ 1 agrees with 
the observed value of r e ~ 0.09 (Koma tsu et al.|20l"T) , in the ab- 
sence of an X-ray contribution to reionization. Although our mod- 
els depend only on the product in eq. (9J, we show on the RHS 
some reasonable values for the component terms. iV 7 = 4400 is 
expected for PopII stars (e.g. Barkana & Loeb 2005}, and studies 
of the high-redshift Lyman alpha forest suggest n rec ~ ( Miralda-| 



|Escude|2003l|Bolton & Haehnelt|2007l |McQuinn et al.|2011) .On 

the other hand, the parameters / esc and /* are extremely uncer- 
tain in high-redshift galaxies (e.g. |Gnedin et a l. 2008 ; Wise & Cen| 
2009 1 |Ferrara & Loeb|20T2]>, though their product m ight be more 
robust IWise & Cen|2009||Paardekooper et al.|2011| >. 

Similarly, the X-ray efficiency is parameterized as follows: 



fx 



Nx 
0.4 



o.i 



(10) 



The value of Nx = 0.4 X-ray photons per stellar baryon cor- 
responds to a total X-ray luminosity above hvo = 0.3 keV of 
£x,o.3+keV ~ 10 40 erg s _1 (M© yr -1 ) -1 , using our spectral en- 
ergy index of a = 1.5. This choice is consistent with (a factor 
of 1.8 higher than) an ex trapolation from the 0.5-8 keV measure- 
ment of ^^eoe^j2012k L x ,o. 5- skev « 3 x 10 39 erg s -1 (M© 
yr~ 1 )~r ]Mineo et al.| ( |2012| find a factor of 2-3 intrinsic scatter in 
this relation for their sample of star-forming galaxies, while other 
studies in different bands have estimates which are also within a 
factor of few of this valu e (e.g. |Gilfanov et al.|200 4 Ranalli et al. 



2003 Per sic et al.|[2004| [Persi c & Rep haeli||2007| [Lehmer et~ 
2010 [I] Recently, the 4 Ms exposure of the Chandra Deep Field- 
South (CDF-S; e.g. |Xue etal.|2011| ) has allowed us to study the red- 
shift evolution of this relation. There is some evidence of a (weak) 
evolution (Basu-Zych et al. 2012, though see |Cowie et al.|20"T2] >. 
However, interpreting the derived upper limits on the X-ray flux of 
z ~ 6 galaxies is very uncertain, and in particular is heavily de- 
pendent on the assumed extrapolation to the soft X-rays which are 
the ones actually interacting with the IGM. For example, the 2a 
upper limits from |Cowie et al. (2012) still allow Nx ~ 1000 for 
their z ~ 6 galaxy sample, e.g. assuming a SFR of 0.1 M© yr -1 , 
hv = 0.3keVanda = 1.5. 

The strongest constraint on the X-ray contribution to reioniza- 
tion comes from the unresolved soft XRB observed in the CDFs. As 
mentioned in the introduction, the XRB is an indirect probe since 
the observed photons are not the ones interacting with the IGM. 



Therefore, constraints are sensitive to the assumed spectral index or 
any spectral break at high energies. Hickox & Markevitch|{ [2007] ) 



estimate the 1-2 keV unresolved XRB to be 3.4 ± 1.4 x 10 

-2 



-T3~ 



cm 



erg 

deg -2 . This can be compared with the XRB at z = in 
our models, estimated from eq. {6|. X-rays in our fiducial model be- 
low make a negligible contribution to the XRB, contributing 0.4% 
(5%) of the unresolved XRB if the sources remain on until z ~ 10 
(6). On the other hand, in our fx = 10 4 "extreme" model where 
X-rays dominate reionization (see below), the sources at z > 10 
already saturate the unresolved XRB. This is roughly consistent 
with recent simple estimates (c.f. Fig. 1 in McQuinn 2012, see also 

|Dijkstra et al.|2004[|Salvaterra et al.|2007). 

Following previous works (e.g. |Furlanetto|2006[|Pritchard & 



Furlanetto 2007 , Santos et al.|2008] |Warszawski et ai.|2Q09j Baek 



et al. 2010, Mesinger et al. 20lTJ, we assume a simple obscuration 



model for our X-ray sources, with a step function at hvo. This al- 
lows us to separately treat the sources of X-ray and UV photons, 
isolating the imprint of the former. On the other hand, the stellar- 
driven UV is assumed to be soft enough that each photon above the 
Lyman limit produces a single ionization (eq.[8j. Although sources 
of X-rays in principle do also contribute some UV photons, lo- 
cal observations find that the majority of both AGN and HMXBs 
are significantly obscured (e.g. Lutovinov et al. 2005 , Tozzi et al. 
2006 ). Photons below our fiducial choice of hvo = 300 eV have 
optical depths exceeding unity for Nm > 10 21 ' 5 cm -2 , consistent 
with the column densities seen in high-redshift gamma-ray bursts 
(GRBs) ftbtani et al.|2006||Greiner et al.|2009) . 

The final free parameter we vary is T v i r , the minimum virial 
temperature of halos hosting galaxies. Although our formalism al- 
lows for disparate values of T v - ir for X-ray and UV sources, for 
simplicity we use the same value for both populations. Our fidu- 
cial choice of T v i r = 10 4 K corresponds to the atomic cooling 
threshold. Although the first stars are likely hosted by smaller ha- 



Haiman et al. 1996, Abel et al. 



los with T vir ~fewxl(r K (e.g. 
2002 , Bromm et al. 2002), star formation inside such small halos 
was likely inefficient (with a handful of stars per halo), and was 
likely suppressed by heating and feedback proce sses ([Haiman et al. 



2000 , IRicotti et al.|200H|Mesinger et al.|2006l|Haiman & B ryan 
|200^[TseTiakhovich & HiratapOlOHO'Leary & McQu inn|20T2] L 
We therefore do not include radiation from such so-called miniha- 
los. T v i r could also have been larger than 10 4 K: radiative and/or 
mechanical feedback (e.g. Springel & Hernquist 2003) eventually 
suppresses star formation inside T v i r < 10 b K halos, though the 
details and timing of these processes are not well understood at 
highredshifts (e.g |Mesinger & Dijkstra|2008||Okamoto et al.|2008| 
Pawlik & Schaye|200 9 f] 



In summary, the total (integrated) number of photons per stel- 
lar baryon escaping into the IGM in our models is given by A^ 7 / eS c 
for the ionizing UV photons, and iVx for the X-rays. For the UV 
we don't specifically adopt a spectral energy distribution (SED), 
instead assuming that the spectra are soft enough that the majority 
of photons will ionize at most one atom. For the X-rays, we as- 
sume a rest-frame luminosity, L e oc v 1 " 5 , extending from v = uo 
to v — oo. 



6 We note that our fiducial choice (fx = 1) for the X-ray efficiency, Ax, 
is a factor of « 4 lower than that used in prior studies, based on an extrapo- 
lation of the earlier 2-10keV measurement by Gilfanov et al. ( 2004) down 
to huo = 0.2 keV. Therefore our / x = 5 model below is a closer match to 
most prior theoretical studies (e.g. Furlanetto 2006 Pritchard & Furlanetto 
|2007]|Santos et al.|2008||Warszawski et al.|2009||Mack & Wyithe|2011) 



7 It is unlikely that T v i r is much greater than fewxlO 5 K, since these 
values approximately correspond to the faint end of the observed galaxy lu- 
minosity functions at z ~ 6 (e.g.|Bouwens et al.|2 008 Labbe et al.|2010[ 
|Salvaterra et al.|2011[|Finlator et al.|2011) . Additionally, high T vir models 
have difficulty in latching onto the slow evolution of the observed emissiv- 
ity, as inferred from the Lya forest at 3 < z < 6 (see Fig. 12 in Mesinger 
|etal.|2012) . 
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Figure 1. Slices through ionization fields, showing evolution with comoving distance along the light cone (x-axis). All slices are 1.5 Mpc (1 cell) deep. 
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Run name T vir (10 4 K) / uv fx hu (keV) r e X-ray fraction [Af^ y ] 2 (fiK 2 ) 



Tle4_fuvl_fxl_lkeV 


1 


1 


1 


0.3 


0.092 


3% 


2.3 


Tle4_fuvl_fx5_lkeV 


1 


1 


5 


0.3 


0.094 


8% 


2.3 


Tle4_fuvl_fxl00_lkeV 


1 


1 


100 


0.3 


0.107 


46% 


2.1 


Tle4_fuvl_fxl00_lkeV_feedback 


1 -+ 10 


1 


100 


0.3 


0.069 


46% 


1.3 


Tle4_fuv0.5_fx50_lkeV 


1 


0.5 


50 


0.3 


0.092 


45% 


1.8 


Tle5_fuvl_fxl00_lkeV 


10 


1 


100 


0.3 


0.063 


48% 


1.2 


Tle5_fuvl_fxl0000_3keV 


10 


1 


10 4 


0.9 


0.108 


96% 


0.95 



Table 1. Model parameters and related observational signatures, as described in the text. 



We list all of our runs below (summarized also in Table [TJ: 

• Tle4_fuvl_fxl_lkeV - This is our "fiducial" run, with T v i r = 
10 4 K, /uv = 1, fx = 1, and hvo= 0.3 keV, resulting in a 
spectrally- averaged mean photon energy hv — 0.9 keV. This model 
predicts r e = 0.09^] consistent with WMAP7. In this fiducial 
case, UV photons dominate reionization, with X-rays accounting 
for only w 3% of the total ionizations. 

• Tle4_fuvl_fx5_lkeV - This run differs from the fiducial only 
in having an X-ray efficiency which is a factor of ~5 higher, mak- 
ing it more in line with several previous works (e. g. |Furlanetto 
20061|Pritchard & Furlanetto|20071 [Santos et al.|2008l|Warszaw'ski 
et al.|2009||Mack & Wyithe|2011|[VaIdes et al.|2012) . 

• Tle4_fuvl_fxl00_lkeV - Here we increase the X-ray effi- 
ciency by a factor of 100 from our fiducial choice. In this model, 
X-rays ionize approximately half of the IGlvj^] We use this model 
as a reference scenario in which X-rays are important, and vary 
additional parameters below. 

• Tle4_fuvl_fxl00_lkeV -feedback - This model is the same 
as Tle4_fuvl_fxl00_lkeV but includes a simple, extreme pre- 
scription for thermal feedback, with T v i r increasing to 10 5 K, 
when the mean temperature of the IGM surpasses 10 4 K. Al- 
though thermal (i.e. radiative feedback) was initially suspected 
of being important (e.g. |Shapiro et al.||1994| |Thoul & Weinberg 
1996, Gnedin 2000), later studies concluded that its impact was 

8 Interestingly, the fact that reionization is "inside-out" on large scales 
results in a different r e then would be estimated assuming homogeneous 
reionization. This is due to the fact that r e is proportional to the sightline- 
averaged integral of the ionized fraction times the density: (xi x n), and 
(xi x n) / (x^ x (n). In patchy reionization models, where the high- 
redshift overdensities hosting early galaxies ionize first, this inequality re- 
sults in values of r e that are a few percent higher than estimated assuming 
Xi and n are uncorrected (e.g. r e = 0.092 instead of r e = 0.088 in 
our fiducial model). Understandably, this bias decreases as reionization be- 
comes more homogeneous, such that in our "extreme" model r e has the 
same value, 0.11, computed both ways. The correlation of the ionization 
and density fields can have a < 10% effect on other global signals, such 
as the 21cm brightness temperature during the advanced stages of reioniza- 
tion (Wheeler et al., in preparation). We therefore stress that the correlations 
between various cosmic fields should be taken into account for precise esti- 
mates of even globally-averaged signals. 

9 The relative contribution of UV and X-ray photons to reionization can 
also be estimated from their relative number of ionizations per stellar 
baryon: N 1 f esc for the UV, and NxNi for the X-rays. Here Ni is the 
mean number of ionizations per X-ray photon. Depending on the spec- 
trum and ionized fraction, secondary ionizations can drive Ni ^10-20 
early in reionization, transitioning to Ni = 1 late in reionization when 
all of the primary electron's energy gets deposited as heat. Using Ni « 10 
and the parameter values in our Tle4_fuvl_fxl00_lkeV model, we confirm 
that the number of ionizations by X-rays and UV photons are comparable: 
N x Ni ~ 400 - 7V 7 /eso 



small er at high-redshifts (e.g.|Ricotti et al.|2002l|Kuhlen & Madau| 
|2005l|Mesinger & Dijkstra|2008| |Okamoto e t al. |2008) . Since the 
Tle4_fuvl_fxl00_lkeV_feedback model is extreme in that it in- 
cludes an instantaneous transition to a high T v i r value, it serves 
to bracket the expected effects of feedback. 

• Tle4_fuv0.5_fx50_lkeV - X-rays also contribute about 1/2 of 
the ionizations in this model, but the efficiencies of both the UV 
and X-rays have been lowered by a factor of « 0.52. This is done 
in order to delay reionization so that r e is the same as in our fiducial 
run, facilitating comparison. 

• Tle5_fuvl_fxl00_lkeV - This model has fx = 100 as in 
Tle4_fuvl_fxl00_lkeV, but with T vir = 10 5 K throughout cosmic 
time, corresponding to inefficient star formation in smaller dwarf 
galaxies. 

• Tle5_fuvl_fxl0000_3keV - This is our "extreme" model, in 
which high energy X-rays drive reionization (accounting for over 
96% of all ionizations). Specifically, we take fx = 10 4 , T v - ir = 
10 5 K, and hv = 0.9 keV (requiring N m > 10 23 cm -2 , and 
resulting in a spectrally-averaged photon energy of hv — 2.7 keV). 
As mentioned above, the unresolved XRB is saturated by z > 10 
in this model. 



3 RESULTS 

3.1 Evolution and Timing of Reionization 

In Fig. [T] we show slices through ionization fields, showing evo- 
lution with comoving line-of-sight (LOS) distance along the light 
cone (x-axis). X-rays have a negligible contribution to reionization 
in the fiducial model. The morphology remains unaffected as the X- 
ray efficiency, fx, is increased by a factor of five. However, as the 
X-ray efficiency is increased further, one starts to notice several dif- 
ferences in timing and morphology. As expected, reionization oc- 
curs earlier, has less small-scale structure, with a partially-ionized 
homogeneous "haze" developing in the early stages of reionization. 
We quantify these trends below. 

In the top panel of Fig. [2] we plot the reionization histories 
of our runs. Again, there is little change from the fiducial model 
(black solid curve) as the X-ray efficiency is increased by a factor 
of ^ 5 (dotted red curve). However, increasing the X-ray efficiency 
by a factor of 100 (orange short-dashed curve) causes reionization 
to happen earlier, with the midpoint shifting by Az « 1. 

Increasing the contribution of X-rays to reionization also 
causes reionization to happen more gradually, as can be seen when 
comparing the blue dot-short dashed and black solid curves, which 
are normalized to have the same r e (c.f. the top and fifth pan- 
els in Fig. [TJ. This effect is most seen in the extreme case of 
Tle5_fuvl_fxl0000_3keV, where reionization is dominated by X- 
rays. In this case, dxm / dz is roughly constant after xui < 0.8, and 



© 0000 RAS, MNRAS 000, 000-000 



X-rays in the early Universe 1 



0.6 
0.4 
0.2 



0.1 
I* 0.01 

io- 3 
io- 4 



_ 


1 1 1 1 1 1 1 




_ 


- 






- 




/ / 








/ / 




Tle4_fuvl_fxl_lkeV 


- 


/ / 

II: ' 




Tle4_fuvl_fx5_lkeV ~ 
Tle4_fuvl_fxl00_lkeV - 




if / 





Tle4_fuvl_fxlOO_lkeV J 
-feedback 


'- 1 


if 




Tle4_fuv0.5_fx50_lkeV - 


- 


II: ' 
if 1 

/ , 





Tle5_fuvl_fxl00_lkeV J 


—f- 




Tle5_fuvl_fxl0000_3keV - 




r r | i i i i 


1 1 1 


1 1 | 1 1 1 1 

1 












'*•••>» 






1 llllll 1 1 1 


X'-.. 




^ } 












1 1 1 1 1 1 1 


i , , 


, , 1 , "," , , 



10 



15 



20 



25 



Figure 2. Top panel: Reionization histories of the runs listed in Table 
Bottom panel: The mean ionization fraction, x e , generated only by X-ray 
ionizations (ignoring the UV contribution except to compute the IGM opac- 
ity to X-rays). 




Figure 3. Fraction of i^mter — 3 Mpc regions having a given xm. 
All smoothing is done with a top-hat filter. The runs are denoted by 
the same line styles as in Fig. [2] Tle4_fuvl_fxl_lkeV (black solid), 
Tle4_fuvl_fxl00_lkeV (orange short-dashed), Tle5_fuvl_fxl00_lkeV 
(cyan dot long-dashed),Tle5-fuvlJxl00003keV (magenta short long 
dashed). The panels correspond to xm = 0.8, 0.5, 0.2 (top to bottom). 



the reionization history doesn't show the "knee" feature where the 
slope increases in the middle stages as the collapse fraction grows. 
This is due to the fact that as X-ray driven reionization progresses, 
the increasing abundance of galaxies is countered by the decreas- 
ing efficiency of X-ray ionization: /i on ,j —> as xm — > 0, so 
only primary ionizations from X-ray photons can complete reion- 
ization (e.g. Shull & van Steenberg 1985, Furlanetto & Stoever 
|20T0l|Valdes et al.|2010) . " 

Aside from directly ionizing the IGM, the heat input from 
early X-rays can result in thermal feedback: raising the effective 
Jeans mass required for gas to efficiently accrete onto DM ha- 
los (see the discussion in j ]2.3| ). If this process is important, the 
first generations of galaxies would sterilize themselves, and reion- 
ization would be delayed until more massive DM halos emerge 
(see the fourth panel in Fig. [TJ. The extreme limit of this sce- 
nario is explored in Tle4_fuvl_fxl00_lkeV_feedback (green long 
dashed curve), and results in a reionization history with a tail of 
xm « 90% extending to high-z, as also noted by |Ricotti & Os-| 
P^ ( [2004l ). 

In the bottom panel of Fig. [2] we show the mean ioniza- 
tion fraction generated only by X-ray ionizations (i.e. ignoring the 
contribution from UV photons). In fiducial reionization scenarios 
driven by UV ionizations, x e corresponds to the electron fraction in 
the "neutral" IGM, outside of the almost fully-ionized HII regions 
local to the sources. As xhi decreases, the resulting increase in the 
X-ray "horizon" (see eq. [TJ is the likely cause of the slight steep- 
ening of dxe/dz (note that the filling factor of HII regions, Qhii, 
is still used to estimate the opacity of the IGM when computing 
the X-ray flux |Mesinger et al.|201l) . |Baek et al.|p0T0] > also noted 
that such a redshifting X-ray background can have a non-negligible 
contribution to X-ray driven reionization. This effect is most evi- 
dent in the models where T v i r ~ 10 4 K halos drive reionization; on 
the other hand, more massive sources on the exponential tail of the 
mass function appear too suddenly to show this feature (e.g. |Lidz| 
|et al.|2007| >. 



3.2 Ionization Morphology 

It has been suggested that abundant X-rays in the early Universe 
could have a large impact on reionization morphology (e.g., see 
the discussion in Haiman 2011 ). From Fig.[T] we can see that this 
is the case only in extreme models, with no "soft" (< 1 keV) X- 
rays. In this section we further quantify the impact of X-rays on 
reionization morphology. 

In Fig. [3] we plot the histograms of the neutral fraction in 
^filter = 3 Mpc regions, at different stages of reionization. The 
fiducial run (black solid curve) results in the most inhomogeneous 
ionization field since it has the smallest contribution of X-rays. 

As the X-ray efficiency is increased by a factor of 100 (orange 
short-dashed curve), the ionization field becomes noticeably more 
homogeneous. There is a significant suppression (a factor of «3) in 
the number of fully ionized regions in the early stages (xhi ~ 0.8). 
At a fixed global neutral fraction, xm, including X-rays decreases 
the relative contribution of fully-ionized HII regions (driven mostly 
by UV sources). Furthermore, the histogram becomes noticeably 
more peaked with a sharp reduction in the number of highly-neutral 
regions. Specifically, there are no 3 Mpc regions having xhi > 
0.95, 0.85, 0.75 at x m « 0.8, 0.5, 0.2, respectively. 

As sources are hosted by rarer, more biased halos, the reion- 
ization field becomes slightly more inhomogeneous (compare the 
orange and cyan curves). However, this difference is not as large 
as the one imprinted by X-ray reionization, since the DM halo bias 
increases only by several tens of percent going from T v i r = 10 4 K 
to Tvir = 10 5 K jMcQuinn et al.|2007 1. This result indicates that 
a significant contribution (~tens of percent) of X-rays to reioniza- 
tion results in a more homogeneous reionization morphology which 
cannot be countered by having more biased sources. 

The magenta curve in Fig. [3] corresponds to the "extreme" 
model, Tle5_fuvl_fxl0000_3keV, in which X-rays complete reion- 
ization. The sources in this model also have an SED with a more 
energetic mean photon energy of hv w 3 keV. As could be guessed 
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Figure 4. Power spectra of the ionization fields for the same models and 
epochs shown in Fig. [3] 

from Fig. [T] the reionization morphology is this model is signifi- 
cantly more homogeneous, with the ionization histograms sharply 
peaked around xm. 

In Fig. [4] we plot the ionization field power spectra ( A^ x = 
k 3 /(27v 2 V) (|(5 xx | 2 )fe, with £ xx = xm/xm — 1) for the same mod- 
els and epochs shown in Fig. [3] The same trends discussed above 
can also be seen in this figure. The fiducial model displays the stan- 
dard "knee" feature, which shifts to larger scales imprinting the HII 
region size in the early stages of reionization, and either the matter 
power spectrum or the mean free path of UV photons (set by Ly- 
man limit systems) in the late stages of reionization (e.g. Alvarez & 
|Abel|2012) . This "knee" feature should be one of the fundamental 
observables of the upcoming 21cm interferometry measurements. 

Increasing the X-ray contribution to reionization to ^50% 
(comparing the orange and black curves) results in a roughly scale- 
free suppression of the power spectrum by a factor of < 2. However 
the "knee" feature is still preserved. Further moving the sources to 
more biased halos (comparing the orange and blue curves) results 
in a more peaked bubble distribution, with a suppression of small- 
scale power and a boost in large-scale power. We note that the 
small-scale power is most sensitive to the X-rays, while the large 
scale power is most sensitive to the halos which host the sources 
(e.g. |McQuinn et al.|2007| ). On the other hand, the extreme model 
with X-ray driven reionization results in a very large (1-2 orders of 
magnitude) suppression of power, especially on small-scales. 

These trends are qualitatively evident in Fig. [5] where we show 
slices through the ionization fields at xhi = 0.9. The first panel 
shows just ionizations from X-rays in the Tle4_fuv0.5_fx50_lkeV 
model (turning off the UV sources). The second panel shows the 
just the ionizations resulting from UV photons in the same model 
(turning off the X-ray sources). It is evident that X-rays result in a 
smoother ionization field. However, even without any UV sources, 
the soft X-rays are able to carve out highly ionized regions sur- 
rounding the first, most highly biased sources. 

The third panel shows the total ionization field in the same 
model, resulting from both X-ray and UV sources (note that this 
field is not just the sum of the first two, but instead corresponds 
to an earlier epoch since all panels are chosen at xhi ~ 0.9). We 
see that small ionization structures are most affected by X-rays, 



as also noted in Fig. [4] When compared at the same xhi, the ion- 
izations contributing to the small, late appearing, UV-driven HII 
regions, effectively shift to the uniform, partially ionized, X-ray 
driven "haze". Larger HII regions are less affected, since these cor- 
respond to highly-biased regions hosting the first X-ray sources. In 
these regions the X-rays themselves pre-ionized a large fraction of 
the HI, making it easier for the UV photons to finish the job. 

The fourth panel corresponds to the Tle5_fuvl _fx 10000 _3keV 
model. Unlike the model shown in the first panel, the 
Tle5_fuvl_fxl0000_3keV model contains no soft X-rays with en- 
ergies below 0.9 keV. The higher energy photons in this extreme 
model understandably create a more uniform ionization field. 



3.3 The kinetic Sunyaev-ZePdovich signal from 
inhomogeneous reionization 

One of the most promising near-term probes of reionization is the 
kSZ signal. The kSZ results from the scattering of cosmic mi- 
crowave background (CMB) photons off of inhomogeneities in the 
ionization and vel ocity fields ([Sunyaev & Zeldovich 1980 ; Ostriker 
|& Vishniac|1986l|Vishniac|19871|Ma & Fry|2002| ). The kSZ can be 
decomposed into the post-reionization, so-called Ostriker-Vishniac 
(OV; |Ostriker & VishniacTT986l |Ma~& Fry|2002| ) component, and 
a component from during reionization referred to as the patchy 
kSZ ( |Gruzinov & Hu|| 1 998 [ [Knox et al.|19"98) . Observational ef- 
forts with the Atacama Cosmology Telescope (ACT) and the South 
Pole Telescope (SPT) are poised to detect this signal for the first 
time, with projected 1 /iK 2 -level sensitivity to the dimensionless 
kSZ power spectrum around a multipole of / = 3000, [A^ooo] 2 . 
/ = 3000 is the "sweet-spot", corresponding to scales small enough 
not be dominated by the primary anisotropics and large enough not 
to be dominated by the CIB. As an integral measurement, [A^ooo] 2 
is sensitive to the duration of the epoch during which ionized struc- 
tures were of comparable size (I = 3000 corresponds to « 20 Mpc 
at high-z). 

Indeed, recent SPT measurements place a bound of 
[Azsooo] 2 < 2.8 pK 2 at 95% confidence limit (C.L.), which de- 
grades to [Aj3ooo] 2 < 6 /jK 2 if a significant correlation between 
the tSZ and the CIB is allowed ( [Reichardt et al.|2012| ). |Mesinger| 



et al. 



( 2012) showed that the [A Z30 oo] 2 < 278 fiK 2 constraint is 
inconsistent with all physically-motivated, UV-driven reionization 
scenarios, even with conservatively low estimates of the OV con- 
tribution. Reionization by UV photons is simply too inhomoge- 
neous on « 20 Mpc scales to accommodate this constraint, with 
the patchy kSZ signal contributing at least 1.5 /jK 2 . 

This implies that either: (i) there is a significant correlation be- 
tween the CIB and the tSZ; and/or (ii) the early stages of reioniza- 
tion occurred in a much more homogeneous manner, perhaps driven 
by X-ray photons with large mean free paths. Upcoming combined 
analyses of microwave and Herschel far infrared data will likely 
test the former scenario. Here we explore the latter. 

The patchy kSZ signal can be decreased by: (i) lowering the 
redshift of reionization; (ii) shortening its duration; and/or (iii) de- 
creasing ionization structure on / = 3000 scale^] X-ray driven 
reionization does the latter. 

In Fig. [6] we plot the kSZ angular power spectrum Cf sz = 



10 Note that a measurement of the patchy kSZ power combined with r e 
does not uniquely determine the redshift and duration of reionization, since 
the ionization morphology affects the shape of the kSZ power, in addition 
to its amplitude ( Mesing er et al.|2012) . 
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Figure 5. Slices through ionization fields at xhi = 0.9. The first panel shows just X-ray ionizations in the Tle4_fuv0.5_fx50_lkeV model (i.e. l-x e 
computed with eq.|2j. The second panel corresponds to ionizations from UV photons in the same model, but without any X-ray sources. The third panel is the 
total ionization field in the Tle4_fuv0.5_fx50_lkeV model. The last panel is the ionization field in the Tle5_fuvl_fxl0000_3keV run, where soft X-rays below 
< 0.9 keV are obscured. All slices are 1.5 Mpc deep, and show only the 0.5< xui <1 range for improved contrast. 



T c 2 mb |^T(^)| 2 , where T cm b = 2.73 K and St is the Fourier trans- 
form of 

at 1 r 

S T = — (u) = ct t / dzc(dt/dz)e~ Te(z) n e u- v , (11) 

and u is the LOS unit vector, ctt is the Thomson scattering cross 
section, r e (z) is the Thomson optical depth to redshift z in the di- 
rection u, v(u, z) is the peculiar velocity, and n e (u, z) is the elec- 
tron number density. Following Mesinger et al. (2012), we only 
show the power sourced from z > 5.6, corresponding roughly to 
the lower limit on the end of reionization (McGreer et al. 2011). 
Power spectrum values at I = 3000 from our models are also 
shown in Table 1. The fiducial model has / = 3000 kSZ power 
fromz > 5.6 of 2.3/iK 2 . 

Interestingly, we see that X-rays have only a modest impact on 
the shape of the power spectrum at / > 1000, unless they dominate 



reionization (magenta curve). This is in line with what we saw in 
the previous section: the effects on the morphology at fixed xhi are 
rather small. 

As before, increasing the X-ray efficiency by a factor of « 5 
from the fiducial value has little effect on this observable. When the 
X-ray efficiency is further increased to account for w 1 /2 of ion- 
izations, and the reionization history is matched to the same r e , the 
/ = 3000 patchy kSZ power is decreased by w 0.5/iK 2 (compare 
the solid black and blue dot dashed curves). On the other hand, if 
the UV and X-ray efficiencies are not decreased to match the same 
r e as the fiducial model, the same relative X-ray contribution (~ 
1/2 of ionizations) results in only a w 0.2/iK 2 decrease in patchy 
kSZ power over the fiducial model (compare the solid black and 
orange short dashed curves), on scales smaller than / w 1000 (^60 
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Figure 6. The kSZ power spectrum sourced from z > 5.6. The points 
denote the recent 95% C.L. on the / = 3000 power by Reichardt et al. 
(2012), assuming no tSZ-CIB correlation (bottom point) and allowing for a 
tSZ-CIB correlation (upper point). The points also include a conservatively- 
low contribution of « 2/iK 2 from the z < 5.6 OV signal (e.g. Sha w et ak| 
|2012l|Mesinger et al.|2012) . 

MpcJ^] In this case, the more uniform component resulting from 
the X-rays is counteracted by shifting reionization to earlier epochs. 
An earlier reionization increases the kSZ signal in Fig. [6] through: 
(i) an increase in the mean density of the IGM in the integrand of 
eq. [TT] and (ii) a larger post-reionization (OV) component to the 
signal down to z = 5.6, the redshift cut we use. 

This is illustrated in the top panel of Fig. [7] where we show the 
cumulative fraction of the z > 5.6 kSZ signal at / = 3000 sourced 
at global neutral fractions less than xhi- In the fiducial model, 
^20% of the total z > 5.6 power comes from post-reionization 
(OV) 5.6 < z < 8, and the other 80% comes from the epoch of 
reionization (patchy kSZ). When the X-ray efficiency is increased 
by a factor of 100, these fractions change to 30% (70%) for the 
OV (patchy) kSZ components. However, the shape of this cumu- 
lative density distribution (CDF) is essentially unchanged. On the 
other hand, as we increase T v i r , a larger fraction of the power is 
sourced from the early stages of reionization. As noted in Mesinger 
|et al.|p012) , increasing T v i r results in a more uniform HII bubble 
size, reducing the small-scale structure. In the later stages of reion- 
ization, it is the small-scale structure (the tail of the HII bubble 
distribution), which dominates the power on / = 3000 (~ 20Mpc) 
scales. This is further evidence that using [A^ooo] 2 to place limits 
on the end of reionization is highly model-dependent. 

Understandably, Tle5_fuvl_fxl0000_3keV has the largest im- 
pact on the patchy kSZ signal. In fact, it is the only one of our mod- 
els which is marginally consistent with the recent aggressive SPT 
bounds on [A^ooo] 2 - From the top panel of Fig.[7] we see that 60% 
of the z > 5.6 power in the model just comes from the OV effect 
post-reionization, 5.6 < z < 9. Since this model is already un- 
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Figure 7. The cumulative fraction of the z > 5.6 kSZ signal at I = 3000 
sourced by redshifts less than z (bottom panel) and corresponding global 
neutral fractions less than xm (top panel). We show the same models as in 
Fig.0 



realistically extreme (see the discussion in j ]2.3| ), we conclude that 
X-rays could not by themselves match the recent aggressive bound 
from SPT. There must be a sizable contribution from the tSZ-CIB 
cross-correlation. 



3.4 Cosmological 21cm Signal 

The most promising probe of early X-rays is probably the red- 
shifted 21cm line from IGM neutral hydrogen. Precursor instru- 
ments like EDGES have already ruled out extremely rapid (Az < 
0.06) reionization models from the all-sky signal (Bow man &| 
Rogers 2010). First generatio n interferometers, like the Low Fre- 
quency Array (LOF AR; |Harker et al.||201 oF] Murchison Wide 
Field Array (MWA; |ringay et al.|2 012pl and the Precision Ar- 
ray for Probing the Epoch of Re-ionization (PAPER; Parsons et al. 
|2010|p^] are coming on-line, promising to measure the 21cm power 
spectrum during 7 < z < 11. Second generation instruments, 
such as the Square Kilometer Array (SKA; Mel lema et al.|2012p" 



and the Lunar University Network for Astrophysical Research (LU- 
NAR^] should have the sensitivity to directly image the topology 
of reionization (low signal-to-noise, large-scale images might even 
be possible with first-generation interferometers; Zar oubi et al.| 
2012 ). In addition to reionization, the 21cm signal is also a sensi- 
tive probe of the thermal history of the IGM pre -reionization, when 
early X-rays have the largest imprint (e.g. Furlanetto [2006| ). 

The 21cm signal is usually represented in terms of the offset of 
the 21cm brightness temperature from the CMB temperature, T 7 



along a line of sight (LOS) at observed frequency v (c.f. Furlanetto 
|et al.|2006) : 



1 1 Behavior on larger scales is difficult to interpret due to cosmic variance 
in the integrand of eq. fTT) : n e u • v (Jelic et al. 2010] Jelic et al., in prepa- 
ration). However, these scales are dominated by the primary anisotropies, 
and so are impossible to be used for the kSZ in the near future. 
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Figure 8. Slices through the predicted 21cm signal for our models. The slices show evolution with comoving distance along the light cone (x-axis). For movies 
of some of these models, please see http://homepage.sns.it/mesinger/21cm_Movie.html 
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where Ts is the gas spin temperature, r UQ is the optical depth at 
the 21-cm frequency z^o, £ n i(x, 2) = p/p — 1 is the evolved (Eule- 
rian) density contrast, H(z) is the Hubble parameter, dv r /dr is the 
comoving gradient of the line of sight component of the comoving 
velocity, and all quantities are evaluated at redshift z = vo/v — 1. 

The spin temperature interpolates between the CMB tem- 
perature and the gas kinetic temperature, Tk- Since the observa- 
tion uses the CMB as a backlight, a signal is only obtained if 
Ts — > Tk - This coupling is achieved through either: (i) collisions, 
which are effective in the IGM at high redshifts, z > 50; or (ii) a 
Lyman alpha background [so-called Wouthuy sen-Field (WF) cou- 
pling; Wouthuysen 1952, Field 1958], effective soon after the first 
sources turn on. 

In Fig. [8] we show slices through the 5Tb fields in our mod- 
els. It is immediately obvious that the 21cm signal is a physics-rich 
probe, encoding information on various processes during and be- 
fore reionization. Although the exact timing varies, all of our mod- 
els show the same, fiducial sequence (c.f. Furlanetto 2006, §2.1 in 
|McQuinn & O'LearypoT^ : 

(i) Collisional coupling: The IGM is dense at high redshifts, so 
the spin temperature is uniformly collisionally coupled to the gas 
kinetic temperature, Tk = Ts < T 7 . Following thermal decou- 
pling from the CMB (z < 300), the IGM cools adiabatically as 
T K oc (1 + z) 2 , faster than the CMB, T 7 oc (1 + z). Thus 5T b is 
negative. This epoch, serving as a clean probe of the mater power 
spectrum at z > 100, is not shown in Fig. [8] 

(ii) Collisional decoupling: The IGM becomes less dense as the 
Universe expands. The spin temperature starts to decouple from the 
kinetic temperature, and begins to approach the CMB temperature 
again, Tk < Ts < T 7 . Thus 5Tb starts rising towards zero. De- 
coupling from Tk occurs as a function of the local gas density, with 
underdense regions decoupling first. Eventually (by z ~ 25), all of 
the IGM is decoupled and there is little or no signal. This epoch 
corresponds to the black region in the far right of the panels in Fig. 

in 

(iii) WF coupling (i.e. Lya pumping): The first astrophysical 
sources turn on, and begin coupling Ts and Tk, this time through 
the WF effect. 5Tb becomes more negative, reaching values as low 
as 5T b ~ 100-200 (depending on the offset of the WF and X- 
ray heating epochs). This epoch corresponds to the black— ellow 
transitiorp^in the panels of Fig. [8] 

(iv) X-ray heating: X-rays heat the IGM, with the spin temper- 
ature now coupled to the gas temperature, Tk = Ts. As the gas 
temperature surpasses T 7 , the 21cm signal changes from absorp- 
tion to emission, becoming insensitive to the actual value of Ts 



17 The Lyman a background has two main contributors: (i) X-ray exci- 
tation of HI; and (ii) direct stellar emission of photons between Lya and 
the Lyman limit. In general, we find that either (ii) dominates, or the two 
are comparable. However in our extreme model (see the bottom panel of 
Fig. [8}, the direct excitation by X-rays is clearly driving the WF coupling, 
imprinting a very interesting spatial structure to the epoch. 
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Figure 9. Redshift evolution of the average 21cm brightness temperature 
offset from the CMB. 

(see eq.[l2]). This epoch corresponds to the yellow— )>blue transition 
in the panels of Fig. [8] 

(v) Reionization: the IGM becomes ionized, a process which is 
inside-out on large scales. The signal again approaches zero. This 
epoch corresponds to the blue— ^black transition in the panels of 
Fig-i 

These stages can be seen in more quantitative detail in Fig. [9] 
where we plot the corresponding evolution of the average 5Tb. This 
figure clearly shows the emission peak and absorption trough. The 
position and height/depth of the mean signal depends on the tim- 
ing of the astrophysical milestones discussed abovj^] The photon 
production rates in our models trace the evolution of the collapse 
fraction, / co ii- Hence, if smaller dark matter (DM) halos host the 
first sources, then WF coupling, X-ray heating and reionization will 
occur earlier, and more graduallj^] 

All other things being equal, the stronger the X-ray efficiency 
of the early sources, the shallower the 5Tb absorption trough. This 
is due to two reasons: (i) X-ray heating can begin before WF cou- 
pling is complete with Ts > Tk', (ii) the ratio T 7 /Tk evolves 
roughly as oc (1 + z) -1 prior to heating, so earlier X-ray heating 
results in a smaller 5Tb contrast even if Ts ~ Tk. 

Following X-ray heating, the (1 — T 7 /Ts) term in eq. |T2| 
goes to unity and the IGM is seen in emission against the CMB. If 
there is subsequently a delay before the onset of reionization, then 
only the cosmological terms in eq. |T2| are not unity, and the mean 
signal decreases as 5% oc y/1 + z. Such an epoch would facili- 
tate a clean measurement of the matter power spectrum at lower, 
more accessible redshifts than the Dark Ages. We see in Fig. [9] 
that this would occur if reionization was driven by more massive 



18 Our fiducial model has a deeper 5Tb absorption trough and a lower 
emission peak than some previous estimates based on Furlanetto ( 2006), 
due to our choice of a lower fiducial X-ray efficiency as discussed in §2.3| 
A lower X-ray efficiency results in a wider WF coupling — > X-ray heating 
transition, and a narrower X-ray heating — > reionization transition. 

19 Even more gradual heating can result in some DM annihilation models, 
since the heating is driven by the collapse fraction in > Mq DM halos 
JValdes et al.|2012| . 
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Figure 10. 21cm power spectra from the same models as in Fig. [3] shown 
at x m = 0.8, 0.5, 0.2 (top to bottom). 
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Figure 11. Ionization-density cross power spectrum for the same models 
as in Fig. [3] shown at #hi = 0.8, 0.5, 0.2 (top to bottom). 

sources than those driving X-ray heating, such as would be the case 
with strong thermal feedback (Tle4_fuvl_fxl00_lkeV_feedback; 
see also the fourth panel in Fig. [8}. 

In Fig. [T0| we show the 21cm power spectra from the same 
models and epochs shown in Fig. [4] Some of the same qualitative 
trends as noted in Fig. [4] can also be seen in Fig. [Toj however the 
differences between the models are smaller. This can be understood 
if we decompose the 21cm power spectrum (post-heating) to first- 
order: 



A 2 21 = fl [AL - 2xm Al d + x 2 m A 2 dd ] , 



(13) 



where Tb is the average brightness temperature in regions with 
xhi = 1, and P is the power spectrum of the ionization (X) and 
density (D) fields. As we saw in j ]3.2[ the more-uniform X-ray 
component decreases A^ x on small scales. However, this mild de- 
crease is compensated by a decrease in the cross-correlation term, 



A^d (shown in Fig. [IT}, which contributes negatively to the 21cm 
power. The impact of the cross term is most evident in the early 
stages of reionization, where the fiducial model (black curve) has 
a factor of < 2-3 less 21cm power at k < 0.1 Mpc -1 than 
Tle4_fuvl_fxl00_lkeV (orange curve). 

we focus on the evolution of the k 



In Fig. 



12 



0.1 Mpc -1 

mode, likely to lie in the middle of the decade of /c-modes probed 
by the first generation interferometers (e.g. |Furlanetto et aLp 009a, 
|Lidz et al.|2008||Chapman et al.|2012] ). The evolution of our fidu- 
cial model agrees very well with the numerical simulations of Baek 
|et al.| ( [20T0| ) (see the solid curves of the S6 model in their Fig. 6; 
this model has a similar X-ray efficiency as our fiducial model), but 
with their evolution shifted to later times (with their resolution they 
were unable to account for the large majority of atomically-cooled 
galaxies). 



As discussed in previous works (Pritchard & Furlanetto 2007 



|Baek et al.|2 010 ; Mesinge r et al.|201 1 



the evolution of the k = 0.1 
Mpc -1 large-scale power has three peaks, corresponding to the 
epochs of Lya pumping, X-ray heating and reionization. The ear- 
lier peaks (especially the X-ray heating one) are larger, due to the 
larger contrast available in the 21cm absorption regime. The posi- 
tion and width of these peaks scale with the model parameters. A 
higher T v i r results in delayed and more rapid transitions, tracing the 
growth of the exponential, high-mass tail of the halo mass function. 
Increasing the X-ray efficiency mainly shifts the X-ray heating peak 
to earlier epochs. If X-ray heating occurs early enough (e.g. fx > 
10), it overlaps with the Lya pumping epoch, decreasing the asso- 
ciated peak in power. Our extreme run, Tle5_fuvl_fxl0000_3keV, 
generally has less power on k = 0.1 Mpc -1 scales throughout 
cosmic time, with the X-ray heating and WF coupling epochs over- 
lapping significantly. The hard X-rays driving both epochs create 
radiation fields with little small-scale structure, while fluctuations 
in both Lya and gas temperature source the large-scale power; this 
results in a negative slope at k — 0.1 Mpc -1 during X-ray heating 
(bottom panel). 

Zooming-in on the reionization epoch, in the right panel of 
Fig. [12] we plot the same quantities as a function of xhi- In the 
advanced stages of reionization, our results are in excellent agree- 
ment with those in Lidz et al. (2008), who only considered reioniza- 
tion under the commonly-used assumption of Ts ^> T 7 (see their 
Fig. 2). |Lidz et al.| (2008) note that the inside-out nature of reion- 
ization causes a drop in large-scale 21cm contrast at early times 
(^hi > 0.9), before the growth of ionization structures pushes the 
power back up. Stated differently in terms of eq. |T3] >, the 21cm 
power before reionization (and assuming Ts ^> T 7 ) traces the 
matter power spectrum, A^d- In the early stages when HII regions 
are small, the large-scale power drops due to the increase in the 
cross-correlation term, A^ d , which contributes negatively. Then at 
xm < 0.9, the HII bubbles grow, dominating the 21cm power 
through the A^ x term. The k — 0.1 Mpc -1 21cm power sub- 
0.5 (|Lidz et al.|[2008| |Friedrich et al. 



sequently peaks at xm 



, with weak scale dependence (i.e. roughly zero slope; |Zahn 



et al. 201 1 Mesinger et al. 2011). These trends are evident in the 
right panel of Fig. [12] 

Including X-rays impacts this scenario in two main ways. 
Firstly, the long mean free path of X-rays, combined with their 
lower ionization efficiency in regions with a higher ionized fraction, 
results in the uniform, partially-ionized "haze", seen most clearly 
in the bottom panel of Fig.[T] Therefore, models with a large frac- 
tion of X-ray ionizations have a weaker and slower growth of the 
A^ d and A^ x terms. This results in a more extended initial stage 
of reionization (with the large-scale 21cm power falling as A^ d in- 
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Figure 12. Evolution of the 21cm power (top panels) and its slope (bottom panels) at k = 0.1 Mpc 1 . Evolution is shown against redshift on the left, and 
against xui on the right. 



creases, before it becomes dominated by the ionization structure, 
i.e. Ax X ). However, this effect is only really evident in the ex- 
treme Tle5_fuvl_fxl0000_3keV model, where the k = 0.1 Mpc -1 
power continues to fall until xhi ~ 0.5, roughly the epoch when 

andg. 



\ 2 



A^d (compare Figures 



11 



+ softer X-ray photons drive the usual evolution of A 2 i described 
above, with the difference in power between the models being less 
than a factor of 2, even when X-rays contribute ^1/2 of the ioniza- 
tions. 

Secondly, and more importantly, we see that in realistic X-ray 
models the assumption ofTs ^> T 7 is not valid in the early stages 
of re ionization. Instead, the reionization and X-ray heating epochs 
overlap in virtually all of our models, resulting in a strong, scale- 
free drop in power at xm > 0.8-0.9. This overlap is most notable 
in our fiducial model, which has the weakest X-ray heating. Note 
that the peak of the mean brightness temperature, 5Tb, in this model 
occurs at z ~ 12, when reionization was already well underway, 
x m ~ 0.8. 



4 CONCLUSIONS 

Recent interest has focused on X-rays as a possibly significant con- 
tributor to cosmic reionization (e.g. |Haim an 2011, Mi rabel et al.| 
|2011[ |Fontanot et al.|[20T2] >. X-ray driven reionization would be 
fundamentally different from a UV driven one. In particular, reion- 
ization could be more extended with a more uniform morphology. 
In this work, we model the observational signatures of X-rays in 
the early Universe, attempting to provide an intuitive framework 
for interpreting upcoming observations. 

We use the public code, 21CMFAST, to generate seven differ- 
ent models of reionization and pre-heating. Our simulations are 750 
Mpc on a side (currently unachievable by standard RT approaches 
hoping to capture the contribution from the dominant galaxy pop- 
ulation), and span the redshift range 5.6 < z < 40. We vary as- 
trophysical parameters, looking for general trends and physical in- 
sight. In addition to studying the timing and morphology of reion- 
ization, we also predict the kSZ and 21cm signatures in our models. 

We find that by increasing the X-ray contribution, reioniza- 



tion occurs earlier, has less small-scale structure, and develops a 
partially-ionized homogeneous "haze". However, when compared 
at the same xm, the impact is modest. Namely, X-rays yield fewer 
small HII regions, with the IGM between HII regions being par- 
tially ionized. Larger HII regions are less affected, since these host 
the first, highly-biased X-ray sources. In these regions the X-rays 
themselves pre-ionized a large fraction of the HI, making it eas- 
ier for the UV photons to finish the job. Specifically, we find that 
models with a significant (tens of percent) contribution of X-rays 
to reionization exhibit: (i) a lack of fully neutral regions, and (ii) 
a suppression of small-scale (k > O.lMpc -1 ) ionization power 
by a factor of < 2. These changes in the reionization morphology 
cannot be countered by having more biased sources (i.e. a higher 
T v i r ), and are therefore a robust indicator of an X-ray contribution 
to reionization. 

The suppression of small-scale ionization structure by X-rays 
also results in a smaller kSZ signal at / = 3000. In particular, if X- 
rays contribute w 1/2 of ionizations, the patchy kSZ power is de- 
creased by « 0.5 /jK 2 with respect to a UV-only model having the 
same reionization history. The shape of the 10 3 < / < 10 4 power 
spectrum remains largely unaffected. The kSZ signal is reduced 
strongly enough to be marginally consistent with the recent aggres- 
sive SPT constraint (Reichardt et al. 2012 ), only if X-rays fully 
reionize the Universe. Since this model is already unrealistically 
extreme (see the discussion in §2.3| >, it is highly unlikely that X-ray 
reionization can by itself match this bound, which assumes no tSZ- 
CIB cross-correlation. We therefore conclude that there must be a 
sizable contribution from the tSZ-CIB cross-correlation. 

In general, X-rays have a small impact on the 21cm signal 
during the advanced stages of reionization (xhi < 0.8), when com- 
pared at the same xm- The difference in 21cm power is less than a 
factor of ~2, even when X-rays contribute ^1/2 of the ionizations. 
The impact of X-rays on 21cm power is less than that on morphol- 
ogy, since the suppression of the ionization power is compensated 
by a decrease in the density-ionization cross-correlation (which 
contributes negatively to the 21cm power). On the other hand, if 
harder X-rays complete reionization (our "extreme" model), the 
21cm power could be lower by more than an order of magnitude 
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in the later half of reionization, and there would be a local mini- 
mum in the evolution of large-scale power at xhi ~ 0.5 (see top 
right panel of Fig. [12]). 

The situation is different in the early stages of reionization 
(#hi > 0.8), as well as the pre-reionization epochs. During the 
pre-reionization epochs, X-rays govern the timing and duration of 
IGM heating. The large-scale (k « 0.1 Mpc -1 ) 21cm power shows 
three maxima, corresponding to Lya pumping, X-ray heating and 
reionization, with the X-ray heating one being the highest of the 
three. During the early stages of reionization (xhi > 0.8), X-rays 
can also have a large impact on the 21cm power spectrum. This 
is primarily due to the overlap of the X-ray heating epoch and the 
early stages of reionization. The overlap is strongest in our fiducial 
model, which has the weakest X-ray emissivity. In particular, this 
overlap can lead to > 10 times more power at k ~ O.lMpc -1 
than predicted with the commonly-used simplifying assumption of 
T s > T 7 (e.g. |Mesinger & Furlanetto|2007|pdret al.|2008) . 

On the other hand, if thermal feedback was very efficient, the 
early X-ray emitting galaxies would self- sterilize by increasing the 
Jeans mass in the IGM, and reionization would be driven by more 
massive, later appearing sources (e.g. Ricotti & Ostriker 2004). 
Such a scenario would result in an extended period after the com- 
pletion of X-ray heating and before the beginning of reionization. 
In this interim period, the 21cm line would be seen in emission, and 
would offer a clean probe of the matter power spectrum. This probe 
of cosmology occurs at much lower redshifts (making it easier to 
observe) than the commonly-implored Dark Ages at z > 40 (e.g. 
IFurlanetto et al.|2009b) . ~ 

In closing, we enumerate some quantitative observational sig- 
natures which can be used to constrain the first generations of X-ray 
sources: 

(i) an X-ray contribution to reionization of > tens of percent 
would result in a lack of (almost) fully-neutral (xhi > 0.9) regions 
which are ~ few cMpc in scale, during the last half of reionization; 

(ii) a 21cm power spectrum tracing the matter power spectrum 
in emission can indicate strong thermal feedback from early X- 
rays; 

(iii) if reionization is driven by ^1 keV X-rays, we expect the 
large-scale (k « O.lMpc -1 ) 21cm power to fall until around 
xm ~ 0.5, then rise again until xm ~ 0.2, thereby creating a 
local minimum at xm ~ 0.5 where standard models instead pre- 
dict a local maximum (see the top panel of Fig. [12]). 

(iv) if high-z galaxies are more efficient at producing X-rays 
than local ones, the rise in 21cm power associated with X-ray heat- 
ing would occur earlier, at xhi > 0.9. 



We thank Matthew McQuinn for comments on a draft version 
of this paper. 
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